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, We present a new QCD sum rule with high sensitivity to the continuum regions 

> : 

£T) , of charm and bottom quark pair production. Combining this sum rule with existing 
m ■ 

7— I 1 ones yields very stable results for the MS quark masses, m c (m c ) and rhb(rhb). Com- 

parison of our approach with experimental data allows for a robust theoretical error 

m ■ 

estimate. We have also provided a new evaluation of the lifetime of the r lepton, r r , 
Q^' serving as a strong constraint on a s . 

I. INTRODUCTION 

> ■ 

X ■ 

The determination of the fundamental Standard Model (SM) parameters is important in 



its own right. It provides a test of the SM when results from various sources are compared, 
which can foster our understanding of SM dynamics (such as strong QCD effects). This 
may also lead to hints of new physics beyond the SM, when precise values of the SM pa- 
rameters are compared against the predictions of more fundamental theories. For example, 
gauge couplings do not unify within the SM. This gives extra evidence against simple grand 
unification theories (GUTs) such as SU (5) without supersymmetry, in addition to the non- 
observation of proton decay. On the other hand, gauge couplings seem to unify at a scale 
~2x 10 16 GeV in the minimal supersymmetric standard model, which can be interpreted 
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as a hint for supersymmetry as well as GUTs [1]. Most GUTs [2] also predict the mass ratio 
mb/m T . 

It is generally difficult to obtain reliable information on quark masses. The Particle Data 
Group |3[ lists only ranges for their values, indicating a lack of confidence in methods used to 
evaluate them. Indeed, a s is quite large at the mass scales of the bottom and charm quarks, 
questioning the convergence of perturbative QCD (PQCD). Furthermore, non-perturbative 
effects governed by the scale Aqcd ~ 0.5 GeV could be large, thus potentially compromising 
the validity of perturbative calculations. Two types of conditions are known to improve 
the situation: high energy or inclusiveness. As an example for the former, a s and m& can 
be determined at LEP energies using PQCD. This yields a s {M z ) = 0.1200 ± 0.0028 ^ 
with very little theoretical uncertainty. But b(c) quark effects are small, so that rht>(Mz) = 
2.67 ± 0.50 GeV is not well constrained. 

In a recent work we computed a s from r T , by definition an inclusive quantity and 
known to be quite insensitive to effects from non-perturbative QCD (NPQCD) p|. Likewise, 
we used a set of inclusive QCD sum rules to derive values for mAnic) and rhb{rhb). One of 
these sum rules is new, and its use together with existing ones 0Jl| proves to be a powerful 
tool to constrain the continuum region of quark pair production. This will be particularly 
helpful for the case of the b quark for which precise measurements of R(s) (the inclusive 
hadronic cross section normalized to the leptonic point cross section) or of Rb{s) (exclusive 
cross section for bb pairs) are unavailable. 

In section 2, we determine the heavy quark masses based upon the new sum rule in 
addition to known ones. We point out a puzzling discrepancy between theory and the recent 
BES data. The BES data seem to be lower than theoretical predictions by 30% consistently 
across the moments. In section 3, we compute the r lifetime. In section 4, we summarize 
our results. 



II. SUM RULES AND HEAVY QUARK MASSES 

On the basis of an unsubtracted dispersion relation (UDR) it was shown in Ref. Q| 
that knowledge of m c , nib, and a s is sufficient to compute the charm and bottom quark 
contributions to the QED coupling a(y/i = Mz). Conversely, comparison of this UDR 
with the more traditional approaches using a subtracted dispersion relation (SDR) offers 



information on m c and m&. The resulting equation relates an inclusive integrated cross 
section to a difference of vacuum polarization tensors, viz. 
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Eq. (0) defines a continuous set of sum rules parametrized by t, where the limit t — > 
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1 the first moment of H q {t). Similarly, for each higher moment, A4 n , one 
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We now take the opposite limit in Eq. (|T|). t — » 00, and regularize the divergent expression, 
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Q q and n„ a re the quark charge and the number of active flavors. Using expressions derived 




in Refs. (lj,ll5| and taking the limit t —>■ 00, the sum rule (0) becomes: 
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Here, and are the mass and the electronic partial width of resonance R, and R c q 
denote the continuum regions integrated from M = Mb± for b and M = Mpo for c. The 
regularization (JHJ) together with the scale choices m q = m q {m q ) and a s = a s {m q ) eliminates 
(resums) all logarithmic terms in Eq. Unlike in any of the sum rules (J2J), R q ont appears 
unsuppressed in Eq. (0J) so that m q varies exponentially with the experimental information 
on the resonances. We will use Eq. (0J) to constrain the continuum region and work with 
the following ansatz: 
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where now a s = a s (2M), s' = s + 4(m 2 (2M) - M 2 ), and 



d s (2M) /ll n q 



4M 2 ' 



We will use the form in the second line (applying it to all moments) of Eq. (jSJ) with the 
corresponding change in the regularization in Eq. (j3J). This keeps only the leading logarithms 
resumed but allows for an analytical integration. Eq. (J5J) coincides asymptotically with the 
predictions of PQCD for massless quarks and interpolates smoothly between the vanishing 
phase space at the pseudo-scalar threshold and the strong onset of fermion pair production. 
Unlike when PQCD is applied to R(s) directly and relatively close to the resonance region, 
we minimize the exposure to local quark-hadron duality violations by using QCD inclusively 
and by merely requiring stable results across the moments. No claim is being made about 
the local shape of R q — we only need theoretical information about global averages. It 
should be pointed out that the new sum rule (4) and the choice of ansatz © are logically 
independent ingredients. An explicit ansatz facilitates the discussion, but our results are 
essentially independent of the shape of the continuum. 

We use the narrow resonance data [J, J/^, ^(2S) for the c quark and T(IS'), T(2S), 
T(3S) for the b quark, as the only experimental input. The wider resonances in the con- 
tinuum region are assumed to be accounted for by our ansatz (JHJ) because (i) they decay 
almost exclusively into flavored hadrons; (ii) they interfere with the non-resonating part 
of the continuum rendering a common treatment virtually impossible; (iii) the 5-function 
approximation (which is perfect for the narrow resonances) becomes successively worse; (iv) 
the philosophy of our ansatz supposes that it averages over local cross-section fluctuations; 
and (v) we wish to compare E q . Q directly to experimental data on the charm continuum 
region such as from Beijing [l7j . The narrow resonance contribution to the various moments 
is shown in the second column of Table HI The 3rd column gives the continuum contribution, 
and the 4th column shows the totals to be compared with the theoretical moments in the 
last column, viz. 



The C n are known up to 0(a 2 s ) and taken from Refs. |1J, Q] where they were computed for 
arbitrary renormalization scale \i. It seems appropriate to choose \i = rh q (m q ), eliminating 
all logarithmic terms as there is only one scale in the problem. Indeed, the authors of 
Ref. [13J, who have chosen jj, = 3 (10) GeV for the charm (bottom) quark and then evolved 
to fi = m q , report a variation over the first 5 (7) moments of 122 (312) MeV. (For larger 
moments the a s expansion [19] of the gluon condensate contribution 3] breaks down.) Using 
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TABLE I: Results for the lowest moments, At n , defined in Eq. (J2) for n = (t — > oo) and Eq. (J2J) for 
n > 1. The upper (lower) half of the Table corresponds to the charm (bottom) quark. Each moment 
has been multiplied by 10™GeV 2n (10 2n+1 GeV 2n ). The continuum error is from AAg' c = ±1.47. 
The last column shows the theoretical prediction for m c (m c ) = 1.289 GeV, m^rhb) = 4.207 GeV, 
and a s (Mz) = 0.1211, where the uncertainty is our estimate for the truncation error (see text). 

n 

the same moments 13] but choosing fi = rh q instead, we observe a variation of less than 
27 (16) MeV. This impressive improvement clearly overcompensates for the larger a s . We 
will choose \i = m q in the following. As for the theoretical uncertainty associated with the 
truncation of the perturbative series, we use the method suggested in Ref . . In our case 
this yields the error estimate, 
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(Nc = Ca = 4Ci? = 3) corresponding to ±16a 3 /7r 3 in the C n . Comparing the correspond 
estimate against the exactly known coefficients of the first eight moments up to order o? s 14 , 
lisl shows that with \i = m q , 23 of 24 coefficients are within the estimate, while only 
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TABLE II: The left part shows contributions to the charm moments (x 10 n+1 GeV 2n ) from 2M D o < 
^fs < 4.8 GeV, and the right part from 2M D o < y/s < 3.83 GeV. Following Ref. Q], we computed 
the columns labeled BES by subtracting from the threshold data on R(s) the average, R, below 
threshold. (We applied corrections for the leading s-dependence.) The errors combine the statistical 
and uncorrelated systematic ones of R with those in the continuum region and with the common 
systematics (< 3.5%) of the difference. 



one coefficient would have been underestimated by a factor « 1.437. This seems to be a 
reasonable state of affairs for a la error estimate and corresponds to ±20 MeV for rh r (m r ) 

n 

from A4i, while variation of the renormalization scale [13J] assesses this error to only 1 MeV, 
which is optimistic. We show the estimate J7|) in the last column of Table 

The last two columns of that Table would agree within errors even if we had chosen 
significantly smaller variations in A3 and especially A3 (AA^' = ±1.47 accounts for the er- 
ror introduced by our ansatz and is above and beyond the variations induced by the fit 
parameters) . The reason for our more conservative error is shown in Table |H] It shows 
that Eq. (J3j) with A3 = 0.50 reproduces the n dependence of the moments computed from 
recent data by the BES Collaboration remarkably well. However, our method favors 
A3 ~ A3 ~ 1.97, and thus 30 to 40% larger contributions. Note that the quark parton model 



predicts A3 = 1, while from third order massive QCD corrections one expects A3 > 1 
(in agreement with our results). Table HP also compares the BES data to the ^(35*) con- 
tribution Q] in the narrow width approximation. Even assuming that the ty(3S) resonance 
(M WS) = 3.7699 GeV) saturates the cheeps-section in that region, we observe a direct 
experimental 2a discrepancy between Ref. |12j and r<L 3 m = 0.26 ±0.04 keV [3j]. Thus there 
is a discrepancy between perturbative QCD (our set of sum rules) and the BES data, while 



QCD appears to be consistent with the ty(3S) data within la. The BES data seem to be 
lower than theoretical predictions by 30% consistently across the moments. This constitutes 
a great puzzle which needs to be resolved in the future. We may be able to quote smaller 
errors after this situation has been resolved. Nevertheless, the quark masses can still be 
determined precisely enough through the sum rule approach. 

There is a possible contribution from the gluon condensate lal. It is known up to 
, but its actual value is not well known. Its inclusion lowers the extracted quark 
masses, increases A3, and sharpens the discrepancy with the BES data. We can bound its 
value to ;$ 0.07 GeV 4 by demanding n independent results within the uncertainties. We 
use this bound (with a central value of zero) to account collectively for non-perturbative 
uncertainties. They induce errors of about 29 MeV into m c (m c ) (n = 2) and 2.4 MeV into 
rh b (m b ) (n = 6). 

The parametric uncertainties from a s and the quark masses themselves are correlated in 
a complicated way (i) across the moments, (ii) across the two quark flavors, (iii) between the 
theoretical moments and the continuum contribution, and (iv) with each other. In practice, 



all this is accounted for by performing fits to the moments. Heavy quark radiation 
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2jJ is not resonating and problems associated with singlet contributions 
appear only at 0(a^), so these issues should not introduce further uncertainties into our 
analysis. We will present our final results after discussing the r lifetime. 



III. t LIFETIME 



It was pointed out long ago that the total hadronic decay width of the r lepton can be 



reliably computed in the framework of perturbative QCD 



23j. Employing the ratio 



T(r — > v T + Hadrons) 

which is predicted to be N c = 3 in the lowest order and in the absence of Cabbibo mix- 
ing, the perturbative corrections can be obtained from QCD calculations of the two point 
current- current correlator Il^(g) (i.e., the hadronic contributions to the W boson vacuum 



polarization tensor). Electroweak radiative corrections were calculated in Ref. 24j. Non- 
perturbative effects have also been estimated and found to be small [3], thus inducing little 
uncertainty. Therefore R T with its small experimental uncertainty provides a solid venue to 



determine the strong coupling constant precisely. 

^QCD 

can be expressed as a contour integral 7] along \s\ = m\ in the complex s-plane, 
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where the Adler function has been calculated to the third order in as 
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where K 2 = 1.6398 and K 3 = 6.371, assuming three massless flavors of quarks. 

One way to evaluate i?^ CD is to expand «s(— s) perturbatively in terms of asijn T ) with 
the help of the renormalization group equation (RGE) of as- However, this results in a 
series of as(m T ) with poor convergence. It proves to be expedient to keep the following 
contour integrals in Eq. © 
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which are complicated functions of as but well-behaved if numerically integrated with the 
help of the RGE on the complex plane. In general, 
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which we calculate numerically up to 4-loop order in the (3 function 2^|. 
We have included one-loop electroweak radiative corrections 
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where the large log is resumed by the corresponding RGE 



(phase space) corrections 
an expansion in m^/Am 2 
strangeness is 



We have also included QED 



quark condensate contributions, as well as c quark effects in 
In total, the partial width into hadrons with vanishing net 
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For our analysis, the experimental input is the r lifetime, 

290.96 ± 0.59 fs, (15) 



n t l-B s 



evaluating the partial widths into leptons, T e T + T^, as well as T^ d theoretically. The relative 
fraction of decays with AS = —1, B$ = 0.0286 ± 0.0009 |3j, is based on experimental 
data, since the value for the strange quark mass, rh s (m T ), is not well known, and the 
PQCD expansion, Cqq^, proportional to m 2 s converges poorly and cannot be trusted. Cqq^ 
also multiplies the corresponding d terms in r" d , posing the same but numerically less 
important problem there. We so.ved it, by re.ati„ g C& D to the ratio ^17(^,1') = 



2 

mi — 



0.896 ± 0.034 3] (in which to linear order all universal terms cancel), and find Cq C q 
ml) = m2(0.091 ±0.046). 

We computed the world average (jUj) by combining the direct value, r T = 290.6±1.1 fs jjjj], 
with T T (B e ,B^) = 291.1 ± 0.7 fs derived from the leptonic branching ratios B e = 0.1784(6) 
and B^ = 0.1737(6) 3] taking into account their 1% correlation. The dominant theoretical 
error induced by the unknown coefficient c?3 = ± 77 [20( is itself strongly a s -dependent, is 
recalculated in each call within a fit, and induces an asymmetric a s error. 

Other experimental uncertainties arise from [3] m T = 1.77699(28) GeV, \V u d\ = 
0.97485(46), and Bs- Uncertainties from higher dimensional terms in the operator prod- 
uct expansion, OPE, are taken from Ref. [30( and add up to Ar r (OPE) = ±0.64 fs. We 
assume that an uncertainty of the same size is induced by possible OPE breaking effects j^. 
The unknown five-loop /3-f unction coefficient, (3^ = ± 579 20], contributes mainly to the 
evolution of a s (m T ) to a s (Mz) and less to the A^. The sub-leading errors listed in this 
paragraph amount to ±1.2 fs. We find, a s {m T ) = 0.3561^11 and a s {M z ) = 0.1221±^|, 
in excellent agreement with a.JMz) = 0.1200 ± 0.0028 from Z-decays |4| and most other 
recent evaluations of r T |30l l31|. 



IV. SUMMARY 



We have presented a new QCD sum rule with high sensitivity to the continuum regions 
of charm and bottom quark pair production. Combining this sum rule with existing ones 
yields very stable results for the MS quark masses, m c (m c ) and rhfe(mfe). Comparison of our 
approach with experimental data allows for a robust theoretical error estimate. We have also 



provided a new evaluation of the lifetime of the r lepton, r r , serving as a strong constraint 
on a s . Including r T , and the n = 2 and n = 6 moments for the c and b quark, respectively, 
as constraints in a fit to all data |4] yields, 

a s {M z ) = 0.1211^?, 

m c (m c ) = 1.2891°;°!° GeV, (16) 

m b {m b ) = 4.207l°;°3? GeV. 

These results reduce the error Q| in a(Mz) by 25%. 

Finally, we stress again that there is a discrepancy between the BES data on one hand, 
and QCD and the ^(SS) electronic width on the other. The BES data seem to be lower 
than theoretical predictions by 30% consistently across the moments. This constitutes a 
great puzzle. We may be able to quote smaller errors after this situation has been resolved. 
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